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Abstract 

Laser-induced  breakdown  spectroscopy  (LIBS)  provides  an  alternative  chemical  analytical  technique  that  obviates  the  issues  of  sample 
preparation  and  sample  destruction  common  to  most  laboratory-based  analytical  methods.  This  contribution  explores  the  capability  of  LIBS 
analysis  to  identify  carbonate  and  silicate  minerals  rapidly  and  accurately.  Fifty-two  mineral  samples  (18  carbonates,  9  pyroxenes  and 
pyroxenoids,  6  amphiboles,  8  phyllosilicates,  and  11  feldspars)  were  analyzed  by  LIBS.  Two  composite  broadband  spectra  (averages  of  10  shots 
each)  were  calculated  for  each  sample  to  produce  two  databases  each  containing  the  composite  LIBS  spectra  for  the  same  52  mineral  samples.  By 
using  correlation  coefficients  resulting  from  the  regression  of  the  intensities  of  pairs  of  LIBS  spectra,  all  52  minerals  were  correctly  identified  in 
the  database.  If  the  LIBS  spectra  of  each  sample  were  compared  to  a  database  containing  the  other  51  minerals,  65%  were  identified  as  a  mineral 
of  similar  composition  from  the  same  mineral  family.  The  remaining  minerals  were  misidentified  for  two  reasons:  1)  the  mineral  had  high 
concentrations  of  an  element  not  present  in  the  database;  and  2)  the  mineral  was  identified  as  a  mineral  with  similar  elemental  composition  from  a 
different  family.  For  instance,  the  Ca-Mg  carbonate  dolomite  was  misidentified  as  the  Ca-Mg  silicate  diopside.  This  pilot  study  suggests  that 
LIBS  has  promise  in  mineral  identification  and  in  situ  analysis  of  minerals  that  record  geological  processes. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Minerals  are  the  fundamental  chemical  building  blocks  of 
solid  earth  materials.  Identification  of  the  mineral  assemblage 
present  in  a  rock  provides  a  wealth  of  information  about  the  his¬ 
tory  of  the  rock  such  as  the  pressure,  temperature,  and  chemical 
environment  of  its  formation;  the  pressure  and  temperature  con¬ 
ditions  of  any  post- formation  hydrothermal  alteration  or  recrys¬ 
tallization  events;  evidence  for  exsolution;  and  the  extent  of 
weathering  at  the  Earth’s  surface. 
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More  than  3800  distinct  mineral  species  have  been  identified 
on  Earth  [1].  Because  correct  identification  of  minerals  and 
knowledge  about  their  chemical  composition  is  critical  to  an 
understanding  of  the  genesis  and  history  of  any  particular  rock 
body,  several  methods  are  routinely  used  to  identify  minerals.  The 
physical  properties  of  minerals  (hardness,  color,  streak,  luster, 
direction  and  quality  of  cleavage,  habit,  magnetism,  reaction  with 
hydrochloric  acid)  are  used  in  the  field  to  identify  minerals.  To  be 
successful,  these  tests  require  extensive  knowledge  and  experi¬ 
ence  on  the  part  of  the  investigator;  even  then  they  are  not 
necessarily  definitive.  More  sophisticated  techniques  [1,2]  are 
required  for  extremely  fine-grained  minerals  or  to  distinguish 
between  minerals  with  similar  physical  properties  (e.g.  white 
minerals  of  the  same  or  similar  crystal  habit).  Most  minerals  are 
transparent  to  translucent  in  sections  30  pm  thick  (thin  sections) 
and  have  distinct  optical  properties  when  viewed  under  polarized 
light.  Again,  mineral  identification  by  polarized  light  depends  on 
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Table  1 

Formulas  of  carbonate  minerals  used  in  this  study 


Mineral 

Formula 

Hexagonal  rhombohedral  carbonates 

Calcite  (4)a 

CaC03 

Dolomite 

CaMg(C03)2 

Magnesite 

MgC03 

Rhodochrosite  (2) 

MnC03 

Siderite 

FeC03 

Smithsonite  (3) b 

ZnC03 

Orthorhombic  carbonates 

Aragonite  (2) 

CaC03 

Cerussite 

PbC03 

Strontianite 

SrC03 

Witherite 

BaC03 

Monoclinic  carbonates 

Azurite 

Cu3(C03)2(0H)2 

a  Numbers  in  parentheses  indicate  the  number  of  samples  analyzed  in  this 
study.  Unless  specified,  one  sample  was  analyzed. 
b  Smithsonite  samples  included  blue,  green,  and  purple  varieties. 


the  skill  and  knowledge  of  the  micro scopist.  X-ray  diffraction 
(XRD)  measures  the  distances  between  planes  of  atoms  in  the 
crystal  structure,  permitting  accurate  identification  of  even  very 
fine-grained  minerals.  The  fine-scale  chemical  composition  of 
minerals  is  commonly  determined  using  electron  microprobe 
analysis  (EMPA)  of  polished  thin  sections.  EMPA  technology  has 
expanded  during  recent  years  to  measure  trace  element  composi¬ 
tions.  Major  and  trace  element  concentrations  in  minerals  are  also 
determined  using  Laser-Ablation  Inductively  Coupled  Plasma 
Mass  Spectrometry  (LA-ICPMS).  Each  of  these  techniques  pro¬ 
vides  a  unique  and  useful  data  set.  However,  these  laboratory 
techniques  are  time  consuming  and  destructive. 

Across  many  different  areas  of  the  geosciences,  there  has 
been  a  long-standing  requirement  for  a  real-time  technique  for 
in-field  mineral  identification  and  chemical  analysis.  The  devel¬ 
opment  of  portable  XRF  and  Raman  spectrometry  systems  over 
the  past  decade  has  met  this  need  to  a  limited  degree.  Broadband 
Laser-Induced  Breakdown  Spectroscopy  (LIBS)  has  the  po¬ 
tential  to  measure  all  elements  in  rocks,  minerals,  and  soils  in 
real-time  with  no  sample  preparation.  Field-portable  LIBS  sys¬ 
tems  are  in  development  [3,4]. 

Because  field-portable  LIBS  is  on  the  horizon,  and  because 
there  are  many  potential  mineralogical  applications  for  LIBS,  it 
is  appropriate  to  examine  the  capability  of  this  technique  to 
identify  minerals.  LIBS  is  one  of  the  analytical  technologies 
planned  for  next  generation  of  Mars  rovers  [5];  the  use  of 
standoff  LIBS  will  expand  the  geographical  area  that  can  be 
chemically  analyzed  from  the  rover. 

On  Earth,  many  prospecting  techniques  rely  on  the  iden¬ 
tification  of  certain  key  mineral  assemblages  or  chemical  signa¬ 
tures  in  minerals  as  “pathfinders.”  For  instance,  kimberlitic 
diamonds  are  commonly  associated  with  significantly  more 
abundant  Cr-diopside  pyroxenes  and  garnets  with  high  Cr/Ca 
ratios  [6,7].  Many  rock  units  exhibit  spatial  variations  in  mineral 
composition  that  provide  clues  to  the  processes  that  formed  or 
recrystallized  the  rock.  For  example,  Parker  and  Nicholson  [8] 


used  As  variations  in  geothermal  sinters  (amorphous  silica)  as  a 
guide  to  gold  enrichment.  Currently,  these  variations  are  mapped 
by  collecting  samples  along  one  or  more  traverse(s)  across  the 
unit,  followed  by  sample  preparation  and  analysis  in  the  labo¬ 
ratory.  The  mapping  process  could  be  extended  to  cover  the 
entire  three-dimensional  exposure  of  the  rock  unit,  using  field- 
portable  LIBS,  if  LIBS  spectra  are  sufficiently  sensitive  to  mi¬ 
neral  chemistry. 

This  paper  explores  LIBS  analysis  of  two  important  mineral 
groups:  the  carbonates  and  the  silicates.  Samples  from  the 
teaching  collection  at  New  Mexico  State  University  that  are 
representative  of  the  natural  range  of  composition  of  common 
carbonate  and  silicate  minerals  were  analyzed  by  benchtop  LIBS 
at  the  Army  Research  Laboratory  (ARL).  These  data  were  used 
to  determine  whether  LIBS  1)  can  be  used  to  identify  minerals; 
2)  is  able  to  readily  differentiate  between  members  within  a 
group  of  structurally  similar  minerals;  and  3)  is  sensitive  to 
mineral  composition  or  stoichiometry. 

2.  Theoretical:  mineral  structures  and  compositions 

2.1.  Carbonates 

The  carbonates  comprise  a  class  of  minerals  in  which  a  metal 
ion  is  coordinated  by  the  CO3-  carbonate  molecule.  Carbonate 
mineral  structures  (Table  1)  are  either  hexagonal  rhombohedral, 


Table  2 

Formulas  of  silicate  minerals  used  in  this  study 


Mineral 

Formula 

Pyroxenes  and  pyroxenoids 

Aegirine 

NaFeSi206 

Augite 

(Ca,Na)(Mg,Fe,Al)(Si,Al)206 

Diopside  (2)a 

CaMgSi206 

Orthopyroxene  (2) 

(Mg,Fe)Si03 

Rhodonite 

MnSi03 

Spodumene 

LiAlSi206 

Wollastonite 

CaSi03 

Amphiboles 

Actinolite 

Ca2(Mg,Fe)5Sig022(0H)2 

Anthophyllite 

(Mg,Fe)7Si8022(0H)2 

Glaucophane 

Na2  Mg3Al2Si8022(0H)2 

Hornblende 

(Ca,Na)2„3(Mg,Fe,Al)5Si6(Si,Al)2022(0H)2 

Riebeckite 

Na2Fe5Si8022(0H)2 

Tremolite 

Ca2  Mg5Si8022(0H)2 

Phyllosilicates 

Biotite 

K(Mg,Fe)3AlSi3Oi0(OH)2 

Lepidolite 

K(Li,Al)2_3AlSi3O10(O,OH,F)2 

Muscovite  (2) 

KAl2AlSi3O10(OH)2 

Phlogopite 

KMg3  A1S  i3  0 1  o(OH)2 

Serpentine  (2) 

Mg6Si40io(OH)8 

Talc 

Mg3Si4Oi0(OH)2 

Feldspars 

K-feldspar  (3) 

(K,Na)AlSi308 

Plagioclase  (8) 

NaAlSi308-CaAl2Si208 

Mineral  formulas  from  Klein  [1]. 

a  Numbers  in  parentheses  indicate  the  number  of  samples  analyzed  in  this 
study.  Unless  specified,  one  sample  was  analyzed. 
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Fig.  1 .  Correlation  of  two  LIBS  spectra  against  each  other.  The  intensity  of  each  pixel  for  sample  Plagioclase  1  is  plotted  against  the  intensity  of  the  corresponding  pixel 
for  sample  Plagioclase  2.  The  line  is  the  result  of  a  linear  regression  of  the  two  spectra.  The  linear  trend  of  data  above  the  regression  line  records  a  peak  that  has  higher 
intensity  in  sample  Plagioclase  2  than  in  Plagioclase  1 . 


orthorhombic,  or  monoclinic,  depending  on  whether  the  metal 
ion  is  small  (rhombohedral)  or  large  (orthorhombic),  or  whether 
the  carbonate  is  hydrous  (monoclinic).  This  structural  flexibility 
and  the  nearly  ubiquitous  presence  of  the  carbonate  molecule  at 
and  near  the  Earth’s  surface  result  in  the  widespread  occurrence 
of  the  carbonates  in  rocks  formed  by  sedimentary,  hydrother¬ 
mal,  metamorphic,  and  weathering  processes. 

2.2.  Silicates 

Known  as  the  “rock-forming”  minerals,  silicates  are  the 
dominant  minerals  in  the  Earth’s  crust  and  mantle.  Three  groups 
of  silicates  are  considered  in  this  paper:  the  inosilicates,  phyl- 
losilicates,  and  feldpsars  (Table  2).  The  inosilicates  (pyroxenes, 
pyroxenoids,  and  amphiboles)  have  chains  of  silica  tetrahedra 
parallel  to  the  ocrystallographic  axis.  The  phyllosilicates,  or 
layer  silicates,  are  composed  of  alternating  sheets  of  tetrahedra 
and  sheets  of  octahedra.  The  micas  (biotite,  lepidolite,  musco¬ 
vite,  and  phlogopite)  have  K  in  the  interlayer  site  between 
tetrahedral  sheets.  In  other  phyllosilicates  (talc,  serpentine),  the 


interlayer  site  is  absent  or  contains  atoms  other  than  K.  The 
feldspars,  a  large  and  diverse  group  of  minerals,  are  divided  into 
two  types:  the  K-feldspars  ((K,Na)AlSi308)  and  plagioclase 
feldspar  (NaAlSi308-CaAl2Si208).  Substitution  of  NaSi  for 
CaAl  in  plagioclase  is  common,  and  practically  any  composi¬ 
tion  between  the  end-member  compositions  can  crystallize  on 
Earth.  This  study  analyzed  eight  samples  in  the  plagioclase  solid 
solution  series. 

3.  Experimental 

3.1.  Analytical  methods 

The  ARL  laboratory  broadband  LIBS  system  consists  of  a 
pulsed  laser,  mirror  and  focusing  lens  for  laser  light  delivery, 
adjustable  sample  stage,  a  fiber  optic  for  plasma  light  collection, 
broadband  spectrometer,  and  desktop  computer  for  system 
control  and  data  collection.  A  100  mJ  laser  pulse  of  10  ns 
duration  from  an  actively  Q-switched  Ultra  Big  Sky  Technol¬ 
ogies  Nd-YAG  laser  operated  at  1064  nm  is  focused  by  a 
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50  mm  focal  length  convex  lens  to  generate  a  high-temperature 
microplasma  at  the  sample  surface. 

The  sample  to  be  analyzed  was  placed  on  a  sample  stage  and  the 
height  adjusted  manually  so  that  the  laser  spark  was  observed  to 
occur  on  the  sample  surface.  The  laser  beam  diameter  was  3  mm, 
with  a  beam  divergence  of  1  mr,  so  that  the  laser  spot  size  was  about 
60  pm  and  the  crater  generated  on  the  sample  was  on  the  order  of 
100  pm  diameter.  A  bundle  of  seven  600  pm  core  diameter  optical 
fibers  was  oriented  to  collect  the  light  emission  from  the  plasma. 
The  fibers  transmitted  the  acquired  light  to  an  Ocean  Optics,  Inc., 
Model  2000  broadband  spectrometer,  having  a  resolution  of 
0.1  nm,  to  separate  the  light  from  different  atomic,  ionic,  and 
molecular  constituents  of  the  plasma.  The  seven  optical  fibers  were 
connected  to  seven  separate  spectrometers  and  gated  CCD  camera 
pairs;  each  pair  covers  one  portion  of  the  198-964  nm  spectral 
range.  The  LIBS  spectra  collection  software  and  the  statistical 
LIBS  library  software  were  provided  by  Ocean  Optics  Inc. 

3.2.  Data  processing 

Broadband  LIBS  spectra  from  198  to  964  nm  were  collected 
from  twenty  laser  shots  for  each  of  52  minerals  (18  carbonates,  9 
pyroxenes  and  pyroxenoids,  6  amphiboles,  8  phyllosilicates,  and 
1 1  feldspars).  Each  laser  shot  was  preceded  with  a  cleaning  shot, 
identical  in  energy  with  the  subsequent  analytical  shot.  The 
samples  are  natural  minerals  and  chemically  zoned  on  a  scale  of 
millimeters  to  centimeters.  Thus,  samples  were  moved  several 
mm  after  each  analytical  shot  to  sample  a  new  area  on  the 
mineral  surface  for  the  next  analysis  in  order  to  capture  a  realistic 
average  mineral  composition.  After  the  data  collection,  two 
composite  spectra  (averages  of  10  shots  each)  were  calculated 
for  each  sample  to  produce  two  databases  each  containing  the 
composite  LIBS  spectra  for  the  same  52  mineral  samples. 

The  broadband  LIBS  spectra  were  compared  to  each  other  by 
calculating  the  correlation  coefficient  ( R )  of  a  linear  correlation 
of  the  13,605  pixel  intensities  from  198  to  964  nm  for  each  pair 
of  spectra  (Fig.  1).  Two  very  similar  LIBS  spectra  would  yield  a 
correlation  coefficient  close  to  1;  The  correlation  coefficient 
decreases  as  the  similarity  between  two  spectra  decreases.  This 
method  provides  a  simple  and  efficient  way  to  compare  many 
spectra  and  takes  advantage  of  the  totality  of  information 
contained  in  the  broadband  LIBS  spectra. 

4.  Results  and  discussion 

4.1.  Carbonates 

The  carbonates  yield  highly  variable  LIBS  spectra  (Fig.  2),  in 
part  because  of  the  high  degree  of  ionic  substitutions  in  both  the 
rhombohedral  and  orthorhombic  isostructural  groups,  with  the 
chemical  compositions  of  the  different  carbonate  minerals  clearly 
reflected  in  their  spectra.  For  instance,  smithsonite  (ZnC03)  is  the 
only  carbonate  with  Zn  lines  in  its  LIBS  spectra  (Fig.  2),  the 
405.8  nm  Pb  line  is  present  for  cerussite  (PbC03),  and  witherite 
(BaC03)  displays  strong  Ba  lines.  By  contrast,  Ca  is  an  element 
that  is  widespread  throughout  the  natural  environment  and,  there¬ 
fore,  is  present  in  most  of  the  analyzed  carbonates,  except  in  the 


smithsonsites  and  magnesite  (MgC03).  Cu  is  present,  as  expected, 
in  azurite  (Cu3(C03)2(0H)2)  but  is  also  detected  in  smithsonite. 
This  is  not  surprising,  because  smithsonite  occurs  in  weathered 
hydrothermal  ore  deposits,  which  commonly  also  are  sources  of 
both  Cu  and  Zn. 

The  sensitivity  of  the  LIBS  spectra  to  carbonate  chemistry 
suggests  that  it  should  be  possible  to  discriminate  between 
carbonate  minerals  using  LIBS.  This  is  demonstrated  in  Table  3, 
where  the  names  and  correlation  coefficients  for  the  five  mi¬ 
nerals  with  the  most  similar  LIBS  spectra  are  listed  for  the  52 
minerals  in  this  study.  All  1 8  carbonate  specimens  were  correctly 
identified,  as  defined  by  being  the  mineral  with  the  highest 
correlation  coefficient.  Six  samples  (2  aragonites  and  4  calcites) 
are  polymorphs  of  CaC03.  For  most  of  these  samples,  each  of 
the  five  best  matches  was  a  CaC03  polymorph  (Table  3).  The 
exception  is  calcite  2,  which  has  pyroxenes  (diopside  — 
CaMgSi206,  and  augite  —  (Ca,Na)(Mg,Fe,Al)(Si,Al)206)  and 
the  carbonate  dolomite  —  CaMg(C03)2  as  the  most  similar 
minerals.  This  sample  is  a  black  calcite,  and  correlates  well  to 
diopside,  augite,  and  dolomite  because  of  the  presence  of  Ca, 
Mg,  Fe,  and  relatively  low  Al.  For  similar  reasons,  dolomite 
correlates  with  calcite,  diopside,  and  augite,  and  talc  correlates 
with  the  magnesian  silicate  minerals  talc  and  serpentine. 

Some  carbonate  minerals  do  not  match  the  spectra  of  other 
carbonates  well.  For  instance,  spectra  of  the  three  smithsonite 
samples  (green,  blue,  and  pink  in  color)  correlate  well  to  each 
other,  despite  the  differences  in  color  (7?  =  0.938 1-0.9945). 
However,  the  correlation  coefficient  for  the  next  most  similar 
sample  is  much  lower  (0.5670-0.6948),  suggesting  that  it 
would  be  difficult  to  determine  smithsonite’ s  identity  if  a 
smithsonite  spectra  were  not  in  the  database.  Similar  relation¬ 
ships  are  observed  for  rhodochrosite  (MnC03),  cerussite 
(PbC03),  strontianite  (SrC03),  witherite  (BaC03),  and  azurite 
(Cu3(C03)2(0H)2).  The  two  averaged  spectra  for  these  minerals 
match  well  (and  for  the  second  rhodochrosite  samples),  but  all 
other  minerals  only  offer  poor  matches  (7^<0.  75). 

4.2.  Pyroxenes  and  pyroxenoids 

The  pyroxenes  and  pyroxenoids  are  isostructural  silicate  mineral 
groups  exhibiting  a  broad  range  of  ionic  substitution.  The 
compositional  differences  are  recorded  in  the  LIBS  spectra 
(Fig.  3),  and  all  nine  of  these  minerals  were  correctly  identified 
(Table  3).  The  chemical  associations  observed  in  the  carbonate  data 
are  seen  in  the  pyroxene-pyroxenoid  spectra  as  well.  For  instance, 
the  Na-Fe  pyroxene  aegirine  was  found  to  be  more  similar  to  Na- 
rich  feldspars  than  to  the  other  pyroxenes,  which  tend  to  be  Na- 
poor.  Similarly,  the  LIBS  spectra  of  the  Li-rich  pyroxene 
spodumene  is  similar  to  that  of  the  Li-phyllosilicate  lepidolite 
(R= 0.9050),  but  not  to  any  other  mineral  in  our  database  (7?  <  0.75). 

4.3.  Amphiboles 

The  amphibole  isostructural  group  of  silicate  minerals  is 
capable  of  incorporating  an  extremely  wide  variety  of  elements, 
as  reflected  in  the  LIBS  spectra  (Fig.  4).  All  six  amphibole 
samples  were  correctly  identified  (Table  3).  All  of  the  analyzed 


Table  3 


Results  of  mineral  identification  experiment 


Sample 

Best  match 

Second  best  match 

Third  best  match 

Fourth  best  match 

Fifth  best  match 

Mineral 

CCa 

Mineral 

CC 

Mineral 

CC 

Mineral 

CC 

Mineral 

CC 

Carbonates 

Aragonite  1 

Aragonite  1 

0.9734 

Calcite  3 

0.9370 

Calcite  4 

0.9289 

Calcite  1 

0.9267 

Aragonite  2 

0.8716 

Aragonite  2 

Aragonite  2 

0.9867 

Calcite  3 

0.9802 

Calcite  4 

0.9744 

Calcite  1 

0.9633 

Aragonite  1 

0.9219 

Calcite  1 

Calcite  1 

0.9963 

Calcite  3 

0.9788 

Aragonite  1 

0.9710 

Calcite  4 

0.9237 

Aragonite  2 

0.9063 

Calcite  2 

Calcite  2 

0.9933 

Diopside  2 

0.9482 

Augite 

0.9136 

Diopside  1 

0.9127 

Dolomite 

0.9091 

Calcite  3 

Calcite  3 

0.9959 

Calcite  1 

0.9714 

Aragonite  1 

0.9483 

Calcite  4 

0.9429 

Aragonite  2 

0.9358 

Calcite  4 

Calcite  4 

0.9914 

Aragonite  2 

0.9501 

Calcite  3 

0.9267 

Calcite  1 

0.9070 

Aragonite  1 

0.9007 

Dolomite 

Dolomite 

0.9787 

Diopside  2 

0.9596 

Calcite  2 

0.9379 

Diopside  1 

0.9175 

Augite 

0.9108 

Magnesite 

Magnesite 

0.9928 

Talc 

0.8817 

Serpentine  2 

0.8317 

Siderite 

0.7825 

Orthopyroxene  1 

0.7054 

Rhodochrosite  1 

Rhodochrosite  1 

0.9916 

Rhodochrosite  2 

0.9652 

Azurite 

0.7043 

Rhodonite 

0.6684 

Siderite 

0.5050 

Rhodochrosite  2 

Rhodochrosite  2 

0.9986 

Rhodochrosite  1 

0.9547 

Azurite 

0.6789 

Rhodonite 

0.6027 

Siderite 

0.4476 

Siderite 

Siderite 

0.9890 

Glaucophane 

0.7907 

Magnesite 

0.7790 

Actinolite 

0.7634 

Talc 

0.7601 

Smithsonite — blue 

Smithsonite — blue 

0.9945 

Smithsonite — purple 

0.9818 

Smithsonite — green 

0.9459 

Plagioclase  5  (Ca-Na) 

0.5670 

Wollastonite 

0.5669 

Smithsonite — green 

Smithsonite — green 

0.9919 

Smithsonite — purple 

0.9481 

Smithsonite — blue 

0.9457 

Plagioclase  1  (Na) 

0.6948 

K-feldspar  1 

0.6877 

Smithsonite — purple 

Smithsonite — blue 

0.9827 

Smithsonite — purple 

0.9791 

Smithsonite — green 

0.9381 

Plagioclase  5  (Ca-Na) 

0.6193 

Wollastonite 

0.6179 

Cerussite 

Cerussite 

0.8941 

K-feldspar  1 

0.5474 

K-feldspar  2 

0.5350 

Serpentine  1 

0.5211 

Plagioclase  1  (Na) 

0.5184 

Strontianite 

Strontianite 

0.9906 

Cerussite 

0.4246 

Wollastonite 

0.3954 

Plagioclase  5  (Ca-Na) 

0.3936 

Tremolite 

0.3838 

Witherite 

Witherite 

0.9957 

Cerussite 

0.4336 

Phlogopite 

0.4182 

Muscovite  1 

0.2969 

Smithsonite — green 

0.2929 

Azurite 

Azurite 

0.9881 

Rhodonite 

0.7413 

Rhodochrosite  1 

0.6537 

Plagioclase  5  (Ca-Na) 

0.6473 

Rhodochrosite  2 

0.6342 

Pyroxenes  and  pyroxenoids 

Aegirine 

Aegirine 

0.9979 

K-feldspar  3 

0.9361 

Plagioclase  4  (Na-Ca) 

0.9289 

Plagioclase  7  (Ca-Na) 

0.9153 

Plagioclase  1  (Na) 

0.8998 

Augite 

Augite 

0.9936 

Diopside  1 

0.9924 

Wollastonite 

0.9494 

Diopside  2 

0.9469 

Calcite  2 

0.9355 

Diopside  1 

Diopside  1 

0.9969 

Augite 

0.9944 

Wollastonite 

0.9448 

Actinolite 

0.9239 

Diopside  2 

0.9222 

Diopside  2 

Diopside  2 

0.9972 

Calcite  2 

0.9343 

Dolomite 

0.9342 

Diopside  1 

0.9180 

Augite 

0.9116 

Orthopyroxene  1 

Orthopyroxene  1 

0.9896 

Orthopyroxene  2 

0.9870 

Serpentine  1 

0.9547 

Anthophyllite 

0.9255 

Serpentine  2 

0.9001 

Orthopyroxene  2 

Orthopyroxene  2 

0.9912 

Orthopyroxene  1 

0.9741 

Anthophyllite 

0.9582 

Serpentine  1 

0.9387 

Tremolite 

0.9054 

Rhodonite 

Rhodonite 

0.9947 

Plagioclase  8  (Ca-Na) 

0.8965 

Plagioclase  5  (Ca-Na) 

0.8911 

Tremolite 

0.8878 

Glaucophane 

0.8619 

Spodumene 

Spodumene 

0.9733 

Lepidolite 

0.9050 

Anthophyllite 

0.7404 

Hornblende 

0.7365 

Muscovite  2 

0.6869 

Wollastonite 

Wollastonite 

0.9957 

Diopside  1 

0.9540 

Augite 

0.9506 

Diopside  2 

0.9323 

Plagioclase  5  (Ca-Na) 

0.9285 
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Amphiboles 

Actinolite 

Actinolite 

0.9834 

Glaucophane 

0.9626 

Tremolite 

0.9506 

Augite 

0.9185 

Diopside  1 

0.9141 

Anthophyllite 

Anthophyllite 

0.9829 

Glaucophane 

0.9183 

Tremolite 

0.9160 

Hornblende 

0.9061 

Orthopyroxene  2 

0.8839 

Glaucophane 

Glaucophane 

0.9840 

Tremolite 

0.9528 

Anthophyllite 

0.9449 

Actinolite 

0.9320 

Orthopyroxene  2 

0.9023 

Hornblende 

Hornblende 

0.9788 

Aegirine 

0.9309 

Plagioclase  7  (Ca-Na) 

0.9055 

K-feldspar  3 

0.9048 

Plagioclase  2  (Na) 

0.8746 

Riebeckite 

Riebeckite 

0.9981 

Plagioclase  7  (Ca-Na) 

0.8925 

Aegirine 

0.8337 

K-feldspar  3 

0.8047 

Hornblende 

0.7374 

Tremolite 

Tremolite 

0.9850 

Plagioclase  6  (Ca-Na) 

0.9516 

Plagioclase  8  (Ca-Na) 

0.9445 

Glaucophane 

0.9288 

Plagioclase  5  (Ca-Na) 

0.9265 

Phyllosilicates 

Biotite 

Biotite 

0.9703 

Phlogopite 

0.9336 

Anthophyllite 

0.8625 

Muscovite  1 

0.8624 

Orthopyroxene  1 

0.8326 

Muscovite  1 

Muscovite  1 

0.9904 

Muscovite  2 

0.8834 

Hornblende 

0.8800 

K-feldspar  3 

0.8437 

Anthophyllite 

0.8304 

Muscovite  2 

Muscovite  2 

0.9809 

K-feldspar  1 

0.9581 

K-feldspar  2 

0.9387 

Anthophyllite 

0.8933 

Serpentine  1 

0.8726 

Lepidolite 

Lepidolite 

0.9758 

K-feldspar  1 

0.8801 

Muscovite  2 

0.8755 

Anthophyllite 

0.8631 

K-feldspar  2 

0.8594 

Phlogopite 

Phlogopite 

0.9952 

Biotite 

0.9233 

Serpentine  2 

0.8513 

Talc 

0.8349 

Anthophyllite 

0.8144 

Serpentine  1 

Serpentine  1 

0.9951 

Orthopyroxene  1 

0.9545 

Orthopyroxene  2 

0.9514 

K-feldspar  2 

0.9272 

Anthophyllite 

0.9012 

Serpentine  2 

Serpentine  2 

0.9970 

Talc 

0.9441 

Orthopyroxene  1 

0.9049 

Orthopyroxene  2 

0.8819 

Phlogopite 

0.8686 

Talc 

Talc 

0.9829 

Serpentine  2 

0.9115 

Magnesite 

0.8552 

Phlogopite 

0.7937 

Orthopyroxene  1 

0.7541 

Feldspars 

K-feldspar  1 

K-feldspar  1 

0.9888 

K-feldspar  2 

0.9766 

Muscovite  2 

0.9300 

Plagioclase  1  (Na) 

0.9296 

Plagioclase  4  (Na-Ca) 

0.9243 

K-feldspar  2 

K-feldspar  2 

0.9903 

K-feldspar  1 

0.9681 

Muscovite  2 

0.9406 

Plagioclase  1  (Na) 

0.9055 

Serpentine  1 

0.9049 

K-feldspar  3 

K-feldspar  3 

0.9813 

Plagioclase  7  (Ca-Na) 

0.9455 

Aegirine 

0.9421 

Plagioclase  4  (Na-Ca) 

0.8784 

Riebeckite 

0.8578 

Plagioclase  1  (Na) 

Plagioclase  1  (Na) 

0.9723 

Plagioclase  4  (Na-Ca) 

0.9642 

Aegirine 

0.9637 

Plagioclase  6  (Ca-Na) 

0.9416 

Plagioclase  2  (Na) 

0.9336 

Plagioclase  2  (Na) 

Plagioclase  2  (Na) 

0.9942 

Aegirine 

0.9327 

Plagioclase  1  (Na) 

0.8908 

Plagioclase  6  (Ca-Na) 

0.8874 

Plagioclase  3  (Na-Ca) 

0.8858 

plagioclase  3  (Na-Ca) 

plagioclase  3  (Na-Ca) 

0.9958 

plagioclase  6  (Ca-Na) 

0.9675 

plagioclase  5  (Ca-Na) 

0.9537 

plagioclase  8  (Ca-Na) 

0.9534 

plagioclase  1  (Na) 

0.9472 

plagioclase  4  (Na-Ca) 

plagioclase  4  (Na-Ca) 

0.9829 

plagioclase  1  (Na) 

0.9787 

K-feldspar  1 

0.9381 

plagioclase  6  (Ca-Na) 

0.9325 

K-feldspar  3 

0.9129 

Plagioclase  5  (Ca-Na) 

Plagioclase  5  (Ca-Na) 

0.9909 

Plagioclase  8  (Ca-Na) 

0.9872 

Plagioclase  3  (Na-Ca) 

0.9660 

Plagioclase  6  (Ca-Na) 

0.9523 

Tremolite 

0.9346 

Plagioclase  6  (Ca-Na) 

Plagioclase  6  (Ca-Na) 

0.9854 

Plagioclase  1  (Na) 

0.9725 

Plagioclase  3  (Na-Ca) 

0.9501 

Plagioclase  8  (Ca-Na) 

0.9397 

Plagioclase  4  (Na-Ca) 

0.9304 

Plagioclase  7  (Ca-Na) 

Plagioclase  7  (Ca-Na) 

0.9990 

Aegirine 

0.9265 

K-feldspar  3 

0.8962 

Riebeckite 

0.8917 

Plagioclase  2  (Na) 

0.8650 

Plagioclase  8  (Ca-Na) 

Plagioclase  8  (Ca-Na) 

0.9694 

Plagioclase  6  (Ca-Na) 

0.9514 

Plagioclase  3  (Na-Ca) 

0.9451 

Plagioclase  2  (Na) 

0.9408 

Plagioclase  5  (Ca-Na) 

0.9334 

a  Correlation  coefficient. 
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Fig.  3.  Average  LIBS  spectra  of  representative  pyroxenes  and  pyroxenoids. 


samples  are  rich  in  Mg,  except  riebeckite  (Na2Fe5Si8022(OH)2), 
which  is  clearly  seen  in  the  intensity  of  the  Mg  279.6  nm  line 
(Fig.  4)  as  well  as  in  the  similarity  between  riebeckite  and  the 
Na- feldspars  rather  than  the  other  amphiboles  (Table  3).  The  five 
most  similar  spectra  for  the  amphiboles  were  from  other 
amphiboles  (42%),  from  the  closely  compositionally  related 
pyroxenes  (25%),  and  from  the  feldspars  (33%),  suggesting  that 
this  method  of  comparing  broadband  LIBS  spectra  is  most 
sensitive  to  mineral  elemental  composition  rather  than  stoichio¬ 
metry.  For  instance,  hornblende  ((Ca,Na)2_3(Mg,Fe,Al)5Si6(Si, 
A1)2022(0H)2)  can  be  considered  to  be  actinolite  (Ca2(Mg, 
Fe)5Si8022(OH)2)  with  the  addition  of  Al  and  Na.  However,  the 
composition  of  hornblende  analyzed  in  this  study  is  more  similar 
to  aegirine  (NaFeSi206)  and  the  feldspars  (7?  =  0.8746-0.9309) 
than  to  actinolite,  Ca2(Mg,Fe)5Si8022(OH)2  (7?  =  0.6477).  An 
unexpected  result  of  LIBS  analysis  of  amphiboles  is  the 
ubiquitous  presence  of  Li  (Fig.  4),  even  in  amphiboles  that 
crystallized  in  geologic  environments  that  are  expected  to  be  Li- 
poor  (anthophyllite  in  metamorphosed  ultramafic  metamorphic 


rocks,  for  example).  The  sensitivity  of  LIBS  to  Li  has  great 
promise  for  the  quantitative  analysis  of  Li  in  Earth  minerals  [9] 
as  well  as  for  discriminating  between  various  rock  types  on  Mars 
(Wiens,  personal  communication,  2006). 

4.4.  Phyllosilicates 

Two  types  of  phyllosilicates  were  analyzed  in  this  study:  the 
micas  with  K+  in  the  interlayer  site  between  the  layers  of  tetra¬ 
hedral  and  octahedral  sheets  (biotite,  muscovite,  phlogopite,  and 
lepidolite),  and  those  minerals  with  a  vacant  or  missing  inter¬ 
layer  site  (serpentine,  talc).  These  differences  are  seen  in  the 
LIBS  data;  all  eight  samples  were  correctly  identified  (Table  3). 
The  presence  of  K+  in  the  mica  interlayer  site  is  clearly  seen  in 
the  LIBS  spectra  (Fig.  5).  The  Mg-rich  nature  of  talc,  serpentine, 
phlogopite,  and  biotite  is  reflected  in  the  LIBS  spectra  for  these 
minerals,  as  is  the  Mg-poor  nature  of  muscovite  and  lepidolite. 
As  observed  in  the  amphibole  data,  Li  is  nearly  ubiquitous, 
present  even  in  talc  and  one  serpentine  sample,  which  crystallize 
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Fig.  5.  Average  LIBS  spectra  of  representative  phyllosilicates. 


in  ultramafic  metamorphic  rocks  that  are  widely  regarded  to 
contain  little  Li. 

4.5.  Feldspars 

The  feldspars  crystallize  across  the  entire  range  of  temperatures 
and  pressures  realized  in  igneous,  metamorphic  and  sedimentary 
rocks  and  therefore  comprise  much  of  the  Earth’s  emst.  The 
structure  of  these  minerals  allows  significant  ionic  substitution  that 
records  details  of  the  changing  geochemical  environment  of  their 
formation.  Thus,  the  ability  to  quickly  distinguish  large  numbers 
of  feldspar  samples  on  the  basis  of  their  chemical  compositions 
would  enhance  the  understanding  of  large  bodies  of  rocks  such  as 
metamorphic  belts  and  plutons  and  be  very  helpful  in  economic 
mineral  exploration.  Representative  feldspar  LIBS  spectra  are 
presented  in  Fig.  6;  all  feldspar  samples  were  correctly  identified 
(Table  3).  The  chemistry  of  the  feldspars  is  expressed  in  the 
observed  LIBS  spectra;  K-feldspar  is  the  only  sample  containing  a 
significant  K  peak  and  the  Na  peak  intensity  is  lower  for  the  calcic 
plagioclase  samples  than  for  the  sodic  plagioclase  samples. 


4. 6.  Mineral  identification  using  LIBS 

Although  it  is  encouraging  that  all  of  the  samples  in  this  study 
were  correctly  identified  by  comparing  two  averaged  LIBS  spec¬ 
tra,  an  important  related  experiment  can  be  realized  by  examining 
the  second  most  similar  spectra  in  Table  3.  This  simulates  the  case 
in  which  a  potentially  similar  mineral,  although  not  the  same 
sample,  exists  in  the  data  set.  In  this  case,  the  existing  data  set  of 
52  minerals  is  insufficient  to  portray  the  variety  of  mineral  com¬ 
positions  observed  on  Earth.  However,  the  exercise  reveals  in¬ 
teresting  aspects  of  LIBS  systematics. 

Using  only  the  second  best  matched  minerals  listed  in  Table  3, 
65%  (34  of  52)  of  the  minerals  were  identified  as  a  mineral  from 
the  same  mineral  family  (carbonate,  pyroxene/pyroxenoid,  am- 
phibole,  phyllosilicate,  or  feldspar).  This  percent  changes  slightly 
from  family  to  family,  from  50%  of  the  amphiboles  to  82%  of  the 
feldspars,  depending  on  the  overall  compositional  variability 
present  in  the  family.  There  were  two  reasons  that  minerals  were 
misidentified  as  a  mineral  of  a  different  family.  The  first  reason  is 
that  the  mineral  contains  a  dominant  element  that  is  not  present  in 
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the  other  minerals  in  the  data  base.  These  matches  have  low 
correlation  coefficients  (<0.80);  examples  are  the  carbonates  ce- 
russite  (PbC03),  azurite  (Cu3(C03)2(0H)2),  and  siderite  (FeC03). 
Smithsonite  (ZnC03)  shows  the  same  behavior;  the  correlation 
coefficients  for  other  smithsonite  samples  are  high  (/?= 0.9381- 
0.9818),  but  the  next  most  similar  spectra  have  low  correlation 
coefficients  (R  <0.70).  Identification  of  these  minerals  would  not 
be  possible  without  having  a  representative  sample  in  the  database. 

Other  minerals  that  were  misidentified  by  the  second  match 
have  relatively  high  correlation  coefficients  (R >0.85).  These 
minerals’  spectra  were  matched  to  minerals  with  similar  ele¬ 
mental  composition  but  different  stoichiometry.  For  instance, 
two  Ca-Mg  carbonates  (calcite  2  and  dolomite)  were  matched 
to  the  Ca-Mg  silicate  diopside,  CaMgSi206,  and  the  Mg  car¬ 
bonate  magnesite  (MgC03)  was  matched  to  the  Mg  silicate  talc 
(Mg3Si4Oi0(OH)2)  There  was  a  high  correlation  between  non¬ 
silicate  carbonate  and  silicate  spectra  because  even  the  major  Si 
(288.16  nm)  and  C  (247.9,  283.7  nm)  peaks  have  low  intensity 
in  minerals  (Fig.  2).  Similarly,  the  phyllosilicates  muscovite  2, 
lepidolite,  and  serpentine  1  were  matched  to  other  silicates  with 
similar  elemental  composition  (K- feldspar  1 ,  K- feldspar  1 ,  and 
orthopyroxene,  respectively).  The  match  between  the  Li  mica 
lepidolite,  (K(Li,Al)2_3AlSi3O10(O,OH,F)2),  and  K-feldspar  1, 
((K,Na)AlSi308),  occurs  because  of  the  K  and  Al  peaks.  Ser¬ 
pentine,  Mg6Si4Oi0(OH)8,  and  orthopyroxene,  (Mg,Fe)Si03, 
are  similar  Mg-rich,  Ca-poor  silicates.  The  two  feldspar  samples 
that  were  misidentified  by  their  second  matches  were  both 
matched  to  the  Na  pyroxene  aegirine,  (NaFeSi206),  because  of 
the  intense  Na  peaks.  These  relationships  suggest  that  mineral 
identification  by  LIBS  is  driven  by  elemental  composition,  as 
would  be  expected.  Because  a  mineral  may  be  misidentified  as  a 
member  of  a  different  mineral  family,  it  is  important  to  recog¬ 
nize  that  LIBS  is  yielding  compositional,  and  not  stoichiometric 
information. 

5.  Conclusions 

This  study  demonstrates  that  broadband  LIBS  spectra  can  be 
used  to  correctly  identify  minerals  from  the  carbonate  and  silicate 
families.  By  comparing  averaged  LIBS  spectra  against  a  library  of 
the  same  mineral  spectra,  100%  of  the  52  minerals  in  the  database 
were  correctly  identified.  If  a  mineral’s  spectrum  was  compared  to 
the  other  5 1  minerals,  there  was  a  65%  chance  that  it  would  be 


identified  as  a  mineral  of  the  same  family  (carbonate,  pyroxene/ 
pyroxenoid,  amphibole,  phyllosilicate,  or  feldspar).  However,  all 
minerals  were  correctly  identified  by  their  dominant  elements 
(Mg-rich  minerals,  Ca-rich  minerals,  etc.),  except  for  those  with 
an  element  that  did  not  exist  in  any  other  mineral  in  the  database. 
This  study  provides  support  for  the  continued  development  of 
LIBS,  especially  portable  LIBS,  for  applications  in  the  geological 
sciences. 
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